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Abstract 

Let S he a. set of n points in the unit square [0, 1]^, one of which is the origin. We construct 
n pairwise interior-disjoint axis-aligned empty rectangles such that the lower left corner of each 
rectangle is a point in S, and the rectangles jointly cover at least a positive constant area (about 
0.09). This is a first step towards the solution of a longstanding conjecture that the rectangles 
in such a packing can jointly cover an area of at least 1/2. 
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1 Introduction 

We consider a rectangle packing problem popularized by Peter Winkler [18 1 119 1 [20]. which has been 
open for decades. It is a one-round game between Alice and Bob. First, Alice chooses a finite point 
set S in the unit square U = [0, 1]^ in the plane, including the origin, that is, (0, 0) G S* (Fig. [l|a)). 
Then Bob chooses an axis-parallel rectangle r(s) C [/ for each point s € S such that s is the lower 
left corner of r(s), and the interior of r{s) is disjoint form all other rectangles (Fig. [l[b)). The 
rectangle r{s) is said to be anchored at s, but r{s) contains no point from S in its interior. It is 
conjectured that for any finite set S C U, (0, 0) E S, Bob can choose such rectangles that jointly 
cover at least half of C/ [U [TE[ [T9l I20j . However, it has not even been known whether Bob can 
always cover at least a positive constant area. It is clear that Bob cannot always cover ^ + £ area 
for any fixed e > 0. If Alice chooses 5" to be a set of n equally spaced points along the diagonal 
[(0, 0), (1, 1)], as in Fig. [uc), then the total area of Bob's rectangles is at most 5 + 5^- There 
was no progress on this problem in the last 20 years, even though it appeared several times in the 
literature. 




(aj (b) (cj 

Figure 1: (a) A set 5 of 6 points in a unit square [0, 1]^, including the origin (0, 0). (b) A rectangle packing 
where the lower left corner of each rectangle is a point in S. (c) Ten equally spaced points along the diagonal 
[(0,0), (1, 1)], and a corresponding rectangle packing that covers roughly 1/2 area. 



Outline. In this paper, we present two simple strategies for Bob that cover at least 0.09121 
area. These are the GreedyPacking and the TilePacking algorithms described below. Both 
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algorithms process the points in the same specific order, namely the decreasing order of the sum of 
the two coordinates, with ties broken arbitrarily (hence (0,0) is the last point processed). 

The GreedyPacking algorithm chooses a rectangle of maximal area for each point in S sequen- 
tially, in the above order. 

The TilePacking algorithm constructs a packing of U with staircase-shaped tiles, and then 
chooses a rectangle of maximal area within each tile independently. We next describe how the 
tiling is obtained. Each tile is a staircase-shaped polygon, with a vertical left side, a horizontal 
bottom side, and a descending staircase connecting them. The lower left corner of each tile is a 
point in S. We say that the tile is anchored at that point. The algorithm maintains the invariant 
that the set of unprocessed points are in the interior of a staircase shaped polygon (super-tile), and 
in addition the anchor and possibly other points are on its left and lower sides. Processing a point 
amounts to shooting a horizontal ray to the right and a vertical ray upwards which together isolate 
a new tile anchored at that point, and the new staircase shaped polygon containing the remaining 
points is updated. Since (0, 0) £ S, TilePagking does indeed compute a tiling of the unit square. 

It will be shown shortly (Lemma [l]) that the GreedyPagking algorithm covers at least as 
much area as TilePagking. Hence it suffices to analyze the performance of the latter. The bulk 
of the work is in the analysis of this simple TilePagking algorithm, which involves geometric 
considerations and a charging scheme. 

Related work. Very little is known about anchored rectangle packing. Recently, Christ et al. |8] 
proved that i/ Alice can force that Bob's share to be less than -, then n > 2 . Our result 
indicates that this condition does not materialize for large r, since Bob can always cover at least a 
constant fraction of the area for any n G N. 

Previous results on rectangle packing typically consider optimization problems, and are only 
loosely related to our work. While our focus here is not on the optimization version of the anchored 
rectangle packing problem, in which the total area of the anchored rectangles is to be maximized 
for a given set S of anchors, our algorithms do provide a constant-factor approximation. 

In the classical strip packing problem, n given axis-aligned rectangles should be placed (without 
rotation or overlaps) in a rectangular container of width 1 and minimum height. This problem 
is APX-hard (by a reduction from bin packing). After a series of previous results {e.g., [16l [T7]). 
Harren et al. [TT] recently found a (5/3 + e)-approximation. Jensen and Solis-Oba [13] devised an 
AFPTAS which packs the rectangles into a box of height at most (1 + e)OPT + 1 for every e > 0. 
Bansal et al. [5] gave a 1.69-approximation algorithm for the 3-dimensional version. 

Further related problems are the 2-dimensional knapsack and hin packing problems. Given a 
set of axis-aligned rectangles and a box B, the geometric 2D knapsack problem asks for a subset 
of the rectangles of maximum total area that fit into B. In contrast, the 2D hin packing asks for 
the minimum number of bins congruent to B that can accommodate all rectangles. Jansen and 
Pradel [15] designed a PTAS for the geometric 2D knapsack problem, although it does not admit 
a FPTAS. The weighted version does not admit an AFPTAS and an approximation algorithm 
by Jansen and Zhang [141 guarantees a ratio of 2 + e for every e > 0. For the 2D bin packing, 
Jansen et al. |12| gave a 2-approximation, and Bansal et al. [3j designed a randomized algorithm 
with an asymptotic approximation ratio of about 1.525, improving the previous ratio 1.691 by 
Caprara [HIT]. However, 2D bin packing does not admit an AFPTAS [11[9]. Finally, we mention 
that Bansal et al. |1| gave a PTAS for the rectangle placement problem, which goes back to Erdos 
and Graham jlOj . Here a given set of axis-aligned rectangles should be arranged (without rotation 
or overlaps) so that the area of their bounding box is minimal. 
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2 Constructing a rectangle packing 

In this section we describe the two strategies for Bob and then compare their performance. 

Ordering the points in S. Let S" be a set of n distinct points in the unit square [0, 1]^ such that 
(0,0) G S. Denote by x(s) and y{s), respectively, the x- and y-coordinates of each point s £ S. 
Order the points in S as si,S2, ■■■ ,Sn such that 

x{sj) + y{sj) < x{si) + y{si) 

for I < i < j < n (ties are broken arbitrarily). See Fig. pi Clearly, we have s„ = (0,0). In 
Greedy Packing, Bob chooses rectangles of maximal area for si, . . . , s„ in this order. 

Greed yPacking. For i = l,...,n, choose an axis-aligned rectangle Vi C [0,1]^ of 
maximal area such that the lower left corner of rj is Sj, and r^ is interior-disjoint from 
any Vj, i > j. 

Recall the partial order, called dominance order, among points in the plane. For two points, 
p = (xp, Up) and q = {xg, yg), we say that p ^ q (in words, q dominates p) if 



Xp f^ Xg 



and 



Vp < Vq- 



With this definition, an axis-aligned rectangle with lower left corner c\ and upper right corner C2 
can be written as {p £^ : ci < p < C2}. In particular, any point in r{s) dominates s. 
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Figure 2: (a) Eleven points in [0, 1]^, sorted by decreasing order of the sum of coordinates, (b) Rectangles 
chosen greedily in this order, (c) The tiling of [0,1]^ induced by the dominance order; and the maximal 
rectangles chosen from each tile. Note that rectangle 8 is smaller than the corresponding greedy rectangle. 



We now define interior-disjoint tiles for the set S = {si, . . . , s„} that jointly cover the unit square 
U = [0, 1]^. For i = 1, 2, . . . , n, let tile tj be the set of points in U that dominate Sj, but have not 
been covered by any previous tile tj, j < i. Formally, let 

ti = {p £ [0, 1] : Si <p and Sj -^ p for all j < i}. 

The tiles are disjoint by definition, and they cover U since the origin is in S. Each tile is a 
staircase polygon with axis-aligned sides, bounded by one horizontal side from below, one vertical 
side from the left, and a monotone decreasing curve from the top and from the right. Observe that 
the axis-aligned rectangle spanned by the lower left corner Si and any point p G ij is contained in 
the tile tj. That means that every maximal axis-aligned rectangle in the tile is incident to the lower 
left corner Sj. We can now describe our second strategy for Bob. 
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TilePacking. Compute the tiling U = UiLi *«• ^'^^ i = I, . . . ,n, independently, choose 
an axis-aligned rectangle rj C ti of maximal area. 

By the above observation, the lower left corner of rj is the lower left corner of the tile ti, which is 
Si G S. 

We now show that GreedyPacking always covers a greater or equal area than TilePagking. 

Lemma 1 For each point Sj G S, GreedyPacking chooses a rectangle of greater or equal area 
than TilePagking. 

Proof. The rectangles chosen by the greedy tiling for Sj, j < i, are all disjoint from the tile ti, 
because every point in r{sj) dominates Sj. Hence, Greedy Pagking could choose any maximal 
axis-aligned rectangle from tile ti, but it may choose a larger rectangle (see, for example, r{ss) in 

Fig.[2]). a 

3 Analysis of TilePacking 

Formally, define p as the worst-case performance of TilePacking. In this section, we show that 
TilePagking chooses a set of rectangles of total area p > 0.09121. In our analysis, we will use 
two variables, /3 > 5 and < A < 1. 

Let Tj C ti be the axis-aligned rectangle of maximal area in tile ti, whose lower left corner is Si, 
and let R = {vi : i = 1, . . . , n} he this set of rectangles. If area(rj) > 0.09121 • area(ti) for every i, 
then our proof is complete. However, the ratio area(rj)/area(tj) may be arbitrarily small because 
the tiles can be arbitrary staircase polygons. 

For /3 > 5, we say that tile ti is a fi-tile if area(rj) < -g area(tj). We will give an upper bound 
F[(3, A) on the total area of /3-tiles for every /3 > 5 and < A < 1. This immediately implies that, 
for every /? > 5, the complement of all /3-tiles cover at least 1 — F{f3, A) area, and the total area of 
all rectangles in R is 

1-F(_/3,A) 



n 



/^area(rj) > p> — ^ — . (1) 



«=i 



In Section [3T] we study the properties of individual /3-tiles, and in Section [3^ we prove an upper 
bound on the total area of /3-tiles (for every /3 > 5), which already gives a preliminary bound 
p > 0.07229, using (II|. By integrating over /3, we improve this bound to p > 0.09121 in Section 3.3 
Finally, Section [3.4| shows that TilePagking runs in 0(n log n) time. 



3.1 Properties of a /3-tile 

Let us introduce some notation for describing a single tile ti (Fig. [sFa)). It is bounded from below 
by a horizontal side, denoted a,, and from the left by a vertical side, denoted bi. The width (resp., 
height) of ti is the length of aj (resp., bi), denoted by |aj| (resp., \bi\). We show next that a /3-tile 
ti has a much smaller area than its bounding box (recall that /3-tiles are defined for /3 > 5). 

Lemma 2 Let /3 > 1 and t be a staircase polygon of height h and width w. If the area of every 
axis-aligned rectangle contained in t is less than area(t)//3, tiien 



area(t) < „ , • hw, (2) 



and this bound is the best possible. 
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Proof. Assume, by translating t if necessary, that the lower left corner of t is the origin. Then the 
bounding box of t is [0,'u;] x [0,/i]. Let area(ti) = /3m^, for some u > 0. Then all vertices of t lie 
strictly below the hyperbola arc f{x) = u^ /x, for < x. See Fig. [3[b). The area of the part of 

[0, w] X [0, h] below the curve f{x) = u^ jx is 



n2 + 



u 



dx 



u?-jh 



+ [u^lnxjl^a/^ =u2(i + lnu;-ln('uV^)) = u^(l + ln(/iit;) - Inu^). 



We have shown that area(t) = fiv? < n^(l + ln(/iu;) — Inu^), or /? < l + ln{hw) —Inu'^. Rearranging 
this inequality yields e^~^v? < hw, which implies Q, as required. 

If we approximate the shaded area in Fig. pfb) with a staircase polygon t of height h and width 
w that lies strictly below the hyperbola, then the area of every axis-aligned rectangle contained in 
t is less than u^, and area(t) can be arbitrarily close to -^ • hw. □ 

Sectors and tips. Fix /3 > 5 and consider a /5-tile tj. Decompose ti into rectangular vertical 
sectors by vertical lines passing though its vertices (Fig. p[c)). Each sector is part of some maximal 
axis-aligned rectangle in tj. Hence the area of each sector is less than area(tj)//3, where near equality 
is possible for the leftmost sector. Similarly, we can decompose ti into rectangular horizontal sectors 
by horizontal lines passing through its vertices, and the area of each sector is less than area(tj)//3. 



f{x) = u^/x 




upper 
tip 



main body 

right tip 



S 



(a) 

Figure 3: (a) A /3-tile ti of width 



(b) ^ "" '"" (c) 

and height \bi\. (b) The portion of the rectangle [0, 



w X 



0, h] below 



the hyperbola arc f{x) — u^/x. (c) The decomposition of ti into vertical sectors. The tips of ti are shaded 

Decompose each /3-tile ti into three parts, called right tip, upper tip, and main body, as follows. 
The right tip of ti is cut off from ti by the right-most vertical line that passes through a vertex of ti 
such that the area of the right part is at least area(ti)//3. See Fig. [SJI^c). Similarly, the upper tip of 
ti is cut off from ti by the upper-most horizontal line that passes through a vertex of tj such that 
the area of the part above is at least area(ij)//3. The remaining part, denoted by t'^, is the main 
body of ti. All three parts are staircase polygons. Both the right and the upper tips are unions of 
some sectors of tj. Since the area of each sector is less than area(tj)//3, the area of each tip is at 
least oarea(ti) but less than garea(ti). In particular, since /3 > 5, the right tip of tj is disjoint from 
the upper tip of t,. 

Let a[ and b[, respectively, be the lower and left side of t[. Note that the topmost horizontal 
side and the rightmost vertical side of t[ each contains some point from S, because each contains a 
reflex vertex of the original tile tj. 



Lemma 3 The width of the right tip of ti is at least \ a[ \ 
of the upper tip of ti is at least \h[\, hence 2|6^| < \bi\. 



hence 2|a'| < |oj|. Similarly, the height 
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Proof. By symmetry, it is enough to prove the first claim. Let r' be the maximal axis-aligned 
rectangle in U whose lower side is a[. Since U is a /3-tile, the area of r' is less than oarea(ij). Recall 
that the area of each tip is at least garea(tj), thus area(r') is less than the area of the right tip 
of ti, and so area(r') is less than the area of the bounding box of the right tip of tj. However, the 
height of r' is strictly greater than the height of the right tip (and its bounding box). Therefore, 
the width of r' , which is \a'-\, is less than the width of the right tip of ti. Thus \a'^\ < |aj| — |a^|, or 
2|a^| < \ai\, as required. □ 

Recall that the areas of the right and upper tips of tj are each less than ■garea(tj). Hence 
area(ij) < J±^aicea{t'j^) . Since U is a /3-tile and t'^ C U, the area of every axis-aligned rectangle 
contained in t'^ is less than area(tj)//3 < area(t'j)/(/3 — 4). Applying Lemma [2] for the main body t[ 
yields 

area(t,) < — ^ • area(tj < — ^ • -^^^35- • la^l • \bi\ = -^ ■ |aj • |6J. (3) 

Tall and wide tiles. We distinguish two types of /3-tiles based on the height and width of 
their main body. A /3-tile U is tall if |a'j| < |6^|; and it is wide if |a^| > |6'J. For wide /3-tiles, 
max(|a'J, |6'J) = |a'J, and ([s]) implies 

area(ti) < ^o3F|a-p. (4) 



Similarly, if a /3-tile ti is tall, then area(tj) < 



gM-5 



un'i 



3.2 Upper bound on the total area of /3-tiles 



In this section we give an upper bound -F(/3, A) (see equation (10) further bellow) on the total area 
of all /3-tiles for every /3 > 5 and < A < 1. It is enough to bound the total area of wide /3-tiles 
by ^F(/3,A). By symmetry, the same upper bound holds for the total area of tall /3-tiles. Let 
VF C {1, . . . n} be the set of indices of the wide /3-tiles. 

To begin, for every tile ti we define two adjacent triangles. Let Aj be the isosceles right triangle 
bounded by aj, the line of slope —1 through Sj, and a vertical line through the right endpoint of 
ttj (see Fig. Qa)). Similarly, let Fj be isosceles right triangle adjacent to hi that lies left of tj. The 
two triangles Aj and Fj are not part of the tiling {i, : i = 1 . . . n}, and they may intersect several 
tiles. A key fact is that these two triangles are empty of points from S in their interior, regardless 
whether the tile ti is a /3-tile or not. 

Lemma 4 Fot every i = 1, . . . ,n, the interior of Aj (resp., Ti) is disjoint from S. 

Proof. Suppose, to the contrary, that there is a point Sj in the interior of Aj. Then Sj is processed 
before Sj, and so the left side bj of the tile tj would cut through the horizontal segment Oj, which is 
a contradiction. Similarly, if Sj lies in the interior of Fj, then the lower side aj would cut through 
the left side 6j. □ 

We charge the area of each wide /3-tile ti to the trapezoid Ai C Aj defined below. Let Ai be the 
set of points in Aj that lie vertically below the segment a[ at distance at most A|a^| from it. See 
Fig. W[c). Using Inequality Q, the area of Ai can be bounded from below as follows in terms of 
the area of tf. 

a«a(A) = M±(LZ^ . A|„:| = M?_^ . |„:P > ^i^f^ ■ „ea(,). (5) 
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Figure 4: (a) Tiles ti,. . . ,ti, and the triangles A^ and Ti. (b) The trapezoid ^^ C A^. (c) Every trapezoid 
Ai is contained in the shaded (hexagonal) region. 

The trapezoids Ai, i £ W, are honiothetic copies of each other, and their area depends only on 
\a[\. Note that the triangles Aj (and also the trapezoids Ai) may extend beyond the boundary of 
U = [0, 1]^ {e.g., in Fig. Qb), Aj extends below U). We show that all trapezoids Ai, i G W, lie in 
a polygon whose area is at most 1 + A(3 — A)/8. 



Lemma 5 Every trapezoid Ai, i G W , lies in a polygon of area 1 + A(3 — A)/8. 

Proof. Every trapezoid Ai lies vertically below a segment a'^, which is part of the lower side of 
tile ti. Therefore, Ai cannot extend beyond the left, right, and upper sides of the unit square U . 
Moreover, we show that Ai is contained in the shaded polygon in Fig. pic). 

Consider a trapezoid Ai, and the corresponding lower side Oj of a /3-tile ti, i £ W. Translate them 
vertically down until aj lies on the x-axis. Then apply a dilation centered at the right endpoint of 
Ui (now on the x-axis), such that the left endpoint of Oj becomes (0,0). Finally, apply a dilation 
centered at (0, 0) such that the right endpoint of Oj becomes (1,0). Observe that through all three 
transformations, the trapezoid Ai remains in the shaded polygon. The shaded polygon is the union 
of U and a trapezoid of area ^(1 + "^-5-) 2 = -^(3 — ^)/^- ^ 



The case of pairwise disjoint trapezoids. If the trapezoids Ai, for all i S W, are pairwise 
disjoint, then we can deduce an upper bound for the total area of wide /3-tiles: By Lemma [5} we 
have 

v^ , , , A(3-A) 8 + 3A- A^ 
Y, area(A,) < 1 + ^ ' = . 



Using Inequalities ^ and ([s]), this implies 



y area(ij 



< 



8 + 3A - A2 



2/3 



3A - A2) /3 



A(2-A)-e/5-5 4A(2-A) e^-^' 



(6) 



The general case of overlapping trapezoids. However, it is possible that the trapezoids Ai, 
i G W , are not disjoint. To take care of overlapping trapezoids, we set up a charging scheme, in 
which we choose a set of "large" pairwise disjoint trapezoids. For every trapezoid Ai, i £ W, denote 
by ii the supporting line of Oj. We say that Ai is above Aj (and Aj is below Ai) if Ai and Aj, i ^ j 
and ii is above ij. We next show that if Ai and Aj overlap, then the trapezoid below the other is 
much larger. 
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Lemma 6 Assume that AiCiAj ^ (!>, for some i,j G W,i ^ j ; and Ai is above Aj . Then ijCiAi C a'- 
and (2- A)|a^| < \a'j\. 

Proof. Refer to Fig. p[a). If Ai n Aj 7^ 0, and line £i is above line ij, then the segment a'- has 
to intersect A^. Note that the left endpoint of a'- is Sj S S, and there is some point from S on 

however, there is no point from S in the interior of Aj. 
Therefore, a'- has to traverse both Ai and Aj, hence ij n Aj C a'-. 

The minimum horizontal cross-section of Aj is (1 — A)|a^|, and the width of the right tip of ij is 
at least \a[\ by Lemma [3J It follows that \a'j\ > \ij n Aj| > (1 - A)|a'J + \a[\ = (2 - A)|a-|. □ 



right tip of ti 




right tip of tj 



(a) 

Figure 5: (a) If Ai and Aj overlap, and Qj lies below a^, then 

intersect Aj from above such that the intervals a'- n A^j^ and a' n A.^^ are disjoint. 




(b) Trapezoids Ai-^ and Ai. 



Charging scheme. We introduce a charging scheme among the trapezoids Ai, i £ W. Initially, 
each trapezoid Ai has a charge of area(j4j). We transfer the charges to a subset of pairwise disjoint 
trapezoids. The transfer of charges is represented by a directed acyclic graph G defined as follows. 
The nodes of G correspond to the trapezoids Ai, i G W. li Ai intersects some other trapezoid 
below, we add a unique outgoing edge from Ai to the trapezoid Aj, j G W, whose top side a'- is 
the highest below a[. Observe that all edges of G are oriented downwards, thus G is acyclic. By 
construction, the out-degree of G is at most one. However, the in-degree of a node in G may be 
higher than one. 



Lemma 7 For every trapezoid Aj, the total area of all trapezoids Ai, i ^ j, with a directed path 



in G to Aj is at most 



1 



(1-A)(2-A) 



■^^a.Ti:ea.{Aj). 



Proof. Fix a trapezoid Aj, j G W. Refer to Fig. Of a). For k > 1, let W^ (^ W he the set of 
indices i of trapezoids Ai that have a directed path of length exactly k to Aj in G. We say that 
the trapezoids Ai, i G W^ are on level k. In particular, the trapezoids Ai, i G W^, on level 1 are 
connected to Aj with a directed edge {Ai,Aj) G G. Let Wj = IJfc>i ^f ^e the set of indices of 
all trapezoids Ai, i ^ j, with a directed path in G to Aj. For each Ai, i G Wj, denote by A* the 
unique trapezoid with {Ai,A'![) £ G. 

By Lemma pi every trapezoid Ai, i £ W^, has width at most |aj| < \a'A/{2 — A), and height at 



most 2^1 1 '^'7 1 • Equality is possible if the lower left corner of Ai coincides with the upper left corner 
Sj of Aj (see Ai^ in Fig. p[a)). The gray triangle in Fig. p[a) is the minimum triangle with base 



Packing anchored rectangles 



a'- that contains the maximal possible trapezoid intersecting Aj from above. By triangle similarity, 
the height of this triangle is t^ | a'- 1 , hence its area is (using Equation ([s] 



-^\a,\ 



A 



2(1 -A) 



I ' |2 



area(^j) 



(1-A)(2-A)- 



(7) 



KI/(2-A) 




(b) 

Figure 6: (a) Ai-^ is the largest possible trapezoid that can intersect Aj from above. All trapezoids with a 
directed path in G to Aj can be translated (without overlaps) into the gray triangle, (b) Disjoint trapezoids 
,Ai^ intersect Aj from above. Each of these trapezoids induces a parallel strip Sj. If the strips Sj^ 

lies in the interior of A^j^ . Dashed lines indicate segments 



Ai^ . . . 
and S, 



intersect, and Aj 



\<- 



is above Ai^, then Si 



We claim that X^ieVK- a'rea(Aj) is at most the area of the gray triangle in Figure pTa). To verify 



the claim, we translate every trapezoid Ai, i G Wj, into the pairwise interior-disjoint trapezoids 
lying in the gray triangle. Each trapezoid will be translated in the same direction, (—1,1), but at 
different distances. In order to control the possible location of the translates, we enclose each Ai, 
i S Wj, in a parallel strip. For every i G Wj, draw lines of slope —1 through the two endpoints of 
a[, and denote by Hj the strip bounded by the two lines. Refer to Fig. p[b). 

First consider the trapezoids Ai, i ^ W^, on level 1. By the definition of G, the trapezoids Ai, 
i G Wj, are pairwise disjoint. Label their strips in increasing order from left to right. We show 



are pairwise interior-disjoint. Suppose to the contrary that Hj^ and 
^2, and such that Ai-^ is above Ai^ (if Ai^ is below Ai^, the strips are 



that the strips Sj, i 

Sj2 intersect, where h < 

obviously disjoint). Since the trapezoids are disjoint, the left endpoint of a'^ 



12' 



G S, lies in the 

interior of the isosceles right triangle bounded by the right side of Ai^^, line ij, and the line of slope 
—1 bounding the strip Hj^ from the right. Since A < 1 and |aj^| — |a^ | > \a[ \ by Lemma pi this 
triangle is contained in Aj^ . However, the triangle Ajj is empty of points from S, and we reached 
a contradiction. 

Applying the above argument for the trapezoids on level k, k = 1,2, . . ., we conclude that all 
trapezoids on level A; -|- 1 are pairwise disjoint, and they induce pairwise disjoint parallel strips. 

We are now ready to describe the translation of the trapezoids Ai, i £ Wj. We translate the 
trapezoids in direction (—1,1) in phases /c = 1,2, .... In phase k, consider each Ai, i £ W^ , 
independently. Translate Ai together with all other trapezoids that have a directed path to Ai by 
the same vector in direction (—1,1) until the lower side of Ai becomes collinear with the upper 
side of A* . Each trapezoid Ai remains in its parallel strip Hj , therefore the trapezoids on the same 
level remain pairwise disjoint. After phase k, there is no overlap between a trapezoid Ai, i £ Wj 
of level k and trapezoids in lower levels. When all phases are complete, all trapezoids are pairwise 
interior-disjoint. 
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It remains to show that after the translation, all trapezoids lie in the gray triangle in Fig. pfa). 
From Lemma p^ and since all translations were done in direction (—1, 1), the trapezoids are on or 
to the right of the line of slope —1 passing through Sj. Also from Lemma l6l the right endpoint of 
a'j is to the right of the right side of every triangle Aj, i G Wh hence to the right of the right side 
of every triangle Aj , i G Wj . For every i G Wj , we have | Oj | — \a[\ < \a'^\ from Lemma 3l and so 

I — I I /I < -T-JT = ^- Hence the upper right corner of every Ai, i ^ W}, is below the line of slope —A 

|"i| |"j| |a^| J 

passing through the right endpoint of a'-. Consequently, the upper right corner of every Ai, i & Wj, 
is below this line. Therefore, after the above translation, all trapezoids Ai, i £ Wj, are contained 
in the gray triangle in Figure ^c) . This verifies the above claim and completes the proof of the 
lemma. □ 

Transfer the charges from all nodes to the sinks in G along directed paths. By Lemma [7| the 
total area charged to a sink Aj is less than 

area(A,) + J^^T^^^^^^^^^^j) = (^ll)(^2-Xf ''''^^^^ ^^^ 

The area of every trapezoid Ai, i £ W, is charged to some sink in G, and the sinks correspond 
to pairwise disjoint trapezoids. We can now adjust Inequality ^ to obtain 

^-^ 3-3A + A2 (8 + 3A-A2) /3 _ (3 _ 3A + a2)(8 + 3A - A^) (3 

,w^ (l-^)(2-^)" 4A(2-A) '^i^~ 4A(1-A)(2-A)2 ^^- ^^^ 

The area of all /3-tiles is less than twice the right hand-side of (pi), namely we can set 



(3 -3A + A2)(8 + 3A- A2) /3 
2A(1-A)(2-A)2 



m^)= '- ::,,':'::j--' --^r m 



Prom M , it follows that the total area of all rectangles in R is 

^1-F{/3,X) 1/ (3-3A + A2)(8 + 3A-A2) /3 
P^ o = o 1 



(3 /3 V 2A(1-A)(2-A)2 e^-^ 

_ 1 (3-3A + A2)(8 + 3A-A2) 1 
~P 2A(1-A)(2-A)2 e^=5- 

Whenever F{f3,X) < 1, this already gives a lower bound of p = ^2(1). We have optimized the 
parameters f3 and A with numerical methods. With the choice of /3 = 12.75 and A = 0.45, we obtain 
an initial lower bound of p > 0.07229. 

3.3 Making /3 a continuous variable 

In this section we further improve the lower bound on the covered area to 0.09121 by making /3 a 
continuous variable and using integration. We define the contribution of each point p G tj C [/, as 
u{p) = area(rj)/area(tj). With this definition, we have 

y area(rj) = / / u{p) dA. 

Let /3o > 5 be a parameter to be optimized later (we will choose /3o = 9.955). Partition the 
interval [/3o,C!o) into subintervals of length e > 0: [/3o,oo) = \JuLi[l3j-i,j3j), where j3j = /3o + js- 
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Denote by Bj C U the union of all /3j-tiles, and let Bj = U\Bj. By definition, we have u(P) > l/Pj 
for every p £ Bj, and so JL^j^. u{p) dA > aiea{Bj)/l3j. 

Observe that the sets Bj form a nested sequence Bq Q Bi Q B2 Q ... ^ U. The total 
contribution of all points in U can be written as 

y area(rj) = / / u{p) dA 
~[ J Jpeu 

= [[ u{p) dA + J2 II <P) d^ 



> 



/pGBo j=i J ■Jp€Bj\Bj.i 

area(i?o) v;^ area(i?j \ -Bj-i) 



/5o ^f^ /3, 

00 

ft ' S ft 



°" — / 1 1 \ 



In Section 3.2, we showed that for any j > 0, area(i?j) < F(/3j, A), where F{f3, A) is given by (10) 



It follows that area(-Bj) > 1 — F{f3j, A), and therefore 



j=0 

l-F{f3o,X) ^ (l-F(/3,-+i,A))-(l-F(/3„A)) 

l-F(/3o,A) ^ F(/3,-,A)-F(/3,-+i,A) 

l-Fi/3o,X) ^ 1 F(/3,+i,A)-F(/3,-,A) 
Letting e go to yields 

1 (3_3A + A2)(8 + 3A-A2) 5/ 1 f^^/l ±l\.^ 

Po 2A(1-A)(2-A)2 V e^o ^ V /3 ' d/3 e/^ ^ ""^ 

1 (3-3A + A^)(8 + 3A-A^) 5/ 1 /- / 1 /?-l\ 

/3o+ 2A(1-A)(2-A)2 V e/^« 4 V/3 ' e/^ ^ ^^ 

^ j_ (3-3A + A^)(8 + 3A-A^) 5 /_ J_ , ^ 1 _ ^ , . 

/3o^ 2A(1-A)(2-A)2 V e/5« 4 e^^ Pe^ 

1 (3-3A + A2)(8 + 3A-A2) 5 f^ 1 



Po 2A(1-A)(2-A)2 J^^ /3e/3 

1 (3-3A + A^)(8 + 3A-A^) 5 
^ 2A(1-A)(2-A)2 "^i(^°)' ^^^) 
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where Ei(x) is the exponential integral 

roc -^ 

For every x > 0, this exponential integral can be approximated by the initial terms of the convergent 
series 

Ei(a;) = —7 — Inx — N^ 



k-k\ ' 
fc=i 

where 7 = 0.57721566 ... is Euler's constant; see [2]. With the choice of /3o = 9.955 and A = 0.452, 
we obtain p > 0.09121. 

Taking into account Lemma [T| we summarize our main result in the following theorem. 

Theorem 8 For any finite point set S C U, (0, 0) G S, the algorithm TilePacking chooses a set 
of rectangles of total area p > 0.09121. Consequently, the same guarantee holds for the algorithm 
GreedyPacking. 

3.4 Runtime analysis 

It is not difficult to show that TilePacking can be implemented in 0(n log n) time and 0{n) 
space in the RAM model of computation. The input is a set S oi n points in the unit square 
U = [0,1]^. Clearly, S can be sorted in 0(n log n) time in non- increasing order of the sum of 
coordinates. Assume that the points are labeled si, S2, . • . , s„ in this order. 

We compute the tiles sequentially in n steps. Let Pi = U \ IJj^i tj be the staircase polygon left 
from U after deleting the first i — 1 tiles. We maintain the x- and y-coordinates of the vertices of 
Pi, respectively, in two binary search trees. In step i, we compute the tile ti by shooting a vertical 
(resp., horizontal) ray from s, until it hits the boundary of Pi. The point hit by an axis-parallel 
ray can be found with a simple binary search in O(logn) time. Once the sides Oj and 5j have been 
determined, we insert the x- and y-coordinates of Si into the search trees, and delete the points 
that are not vertices of Pi subsequently. Each point in S is inserted and deleted at most once, so 
the search trees can be maintained in 0(?ilogn) total time. 

Recall that the generated tiles are staircase polygons, and their reflex vertices are points in S. 
Since every point in S is a reflex vertex in at most one tile, the total complexity of the n tiles 
is 0{n). A tile with k reflex vertices contains exactly k + 1 maximal axis-aligned rectangles, and 
one with the largest area can be selected in 0{k) time. Altogether, we can pick an axis-aligned 
rectangle of maximum area from each of the n tiles in 0{n) total time. Hence TilePacking can 
be implemented in 0(n log n) time. 

4 Conclusion 

We have shown that in the 1-round rectangle packing game, no matter how Alice chooses a flnite 
set of points in S, where (0, 0) G S, Bob can always construct a rectangle packing with rectangles 
anchored at the points in S that cover at least a constant area. Bill Pulleyblank and Peter Winkler 
conjectured that Bob can always cover at least 1/2 of the unit square. While we cannot confirm 
this at the moment, we believe that the performance of our GreedyPacking and TilePacking 
algorithms is significantly better than what we proved here. 

We suspect that the problem of finding the rectangles with maximum total area anchored at the 
given points is NP-hard, but this remains to be shown. Our algorithms certainly achieve a constant 
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approximation ratio 0.09121. No efficient exact algorithm or good approximation was previously 
known. 

Higher-dimensional version. The d-dimensional generalization of the 1-round rectangle packing 
game is also very interesting, and almost nothing is known about it. If 5 is a set of equally spaced 
points along the main diagonal of a d-dimensional unit cube U = [0, l]'^, where (0, . . . , 0) G S, then 
the total volume covered by any anchored d-dimensional axis-parallel rectangle packing is roughly 
1/d. Our GreedyPacking and TilePacking algorithms readily generalize to d-dimensions, but 
the performance analysis does not seem to be easily extendible. 

Multi-round versions. A natural generalization of the problem is the multi-round rectangle 
packing game. One can consider two versions, depending on whether the number of rounds is 
known in advance. In the n-round rectangle packing game, both Alice and Bob know the number 
of rounds. In round i, first Alice places a point Si £ [0, 1]^ somewhere outside of Bob's rectangles, 
and then Bob chooses an axis-aligned rectangle rj C [0, 1]^ with lower left corner at Si and interior- 
disjoint from his previous rectangles. Alice has to choose the origin in one of the n rounds. In the 
unlimited rectangle packing game, the number of rounds (or points) is not known in advance. Each 
round goes exactly as in the n-round version, but the game terminates when Alice decides to put 
a point at the origin and Bob chooses his last rectangle incident to the origin. 

For both versions of the multi-round game, Bob could employ a greedy strategy: for each point 
Si, let Tj be an axis-aligned rectangle of maximum area with lower left corner at Sj that is interior- 
disjoint from all previous rectangles. However, our analysis does not extend to these versions of the 
game. In fact, we can show that the greedy strategy cannot guarantee any constant area for Bob. 
Whether Bob can secure a constant fraction of the area by other means in any of the multi-round 
versions of the game remains open. 

Theorem 9 In both multi-round versions of the rectangle packing game, Bob cannot always cover 
f](l) area with a greedy strategy. 

Proof. We show that for every e > 0, Alice can construct a finite sequence of n points si, . . . , Sn, 
such that Bob can cover at most £ area with a greedy strategy. Essentially, Alice can force Bob to 
choose a rectangle from a |-tile (using at most | of the tile's area) and then fence off the remainder 
of the tile so that it cannot be covered later. Alice can make sure that the total area of these |-tiles 
is arbitrarily close to 1, say 1 — |. Then Bob can cover at most | + (1 — |)-|<e area. We proceed 
with the details. 

We say that a staircase polygon P is a /3-staircase, for some /3 > 1, if the area of every axis-aligned 
rectangle contained in P is at most area(P)//3. By Lemma pi for every /3 > 1, /i > 0, and w > 0, 
one can construct a /3-staircase of height h and width w whose area is roughly f]e^~^hw. 

Alice first computes a packing of the unit square U = [0, 1]^ with — staircases, by successively 
choosing interior-disjoint --staircases of smaller and smaller sizes until their total area is at least 
1 — I (Fig. [7Fa,b)). Specilically, in the current step, given an axis aligned rectangle R of height h 
and width w, a - staircase polygon P with the same height and width is anchored at the lower left 
corner of R, and the remaining space R\P is partitioned into vertical or horizontal sectors to be 
processed (see Section Isl). 

Then she slightly shrinks these staircases, to make them pairwise disjoint, and perturbs them to 
ensure that each --staircase Pi contains a unique axis-aligned rectangle rj of maximal area, and 
Tj has the same width as Pj. The point set S contains, for every — staircase in this packing, all 
vertices of Pi, including the lower left corner Sj. In addition, for every --staircase whose lower left 
corner is not on the left side of U, S also contains a point pi very close to the left side of rectangle 
Tj, as shown in Fig. [TJ^c). 
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Figure 7: (a) An initial --staircase packing, (b) In any empty region (shaded gray), Alice can place an 
additional --staircase, (c) The placement of an extra point pi near the left side of the maximal rectangle 
ri C Pi. 

It remains to determine the order in which Alice reveals the points to Bob. The points associated 
with each --staircase Pi are revealed in a contiguous sequence such that the last two points in each 
sequence are the lower left corner Sj followed by the extra point pi. For the lower left corner Sj, Bob 
has to choose the unique rectangle rj of maximal area in Pi, which is adjacent to the lower side of 
Pi. For the extra point pi, Bob has to choose a tall rectangle of negligible area, which covers the left 
side of Pi. These two rectangles guarantee that no subsequent rectangle can cover any additional 
part of Pi, while area(rj) < | • area(Pj). 

To determine the order of sequences of points associated with the staircases, we define a partial 
order over the staircase polygons. Note that in the initial staircase packing, each tip of every Pi is 
adjacent to the left or lower side of another staircase, or the right or upper side of U. This defines 
a partial order: let Pi -< Pj, if the tip of Pj is adjacent to the left or lower side of Pi; or if the lower 
left corner of Pi is a reflex vertex of Pj. Order the staircases in any linear extension of this partial 
order. This ensures that Bob cannot choose a rectangle intersecting the interior of Pi before Alice 
reveals the lower left corner Sj. □ 

Best versus worst greedy strategy. For a finite set 5 C [0, 1]^, and a permutation (ordering) vr 
of S, we can select anchored rectangles greedily (ties are broken arbitrarily) in the order prescribed 
by vr. One can ask which permutation gives the best or the worst performance for a greedy strategy. 
Our main result, Theorem [sj says that for every n-element point set S, (0,0) G S, we can find in 
0(n log n) time a permutation vr for which the greedy strategy covers 0(1) area. In the worst case. 



greedy covers only o(l) area by Theorem K^ In the best case, however, we will show (Lemma 10) 



that greedy is always optimal for some permutation vr. We say that an anchored rectangle packing 
is Pareto optimal if each rectangle r{s) has maximal area assuming that all other rectangles are 
fixed. It is clear that every optimal solution is Pareto optimal. 

In particular, each rectangle in an optimal solution is bounded by the two rays (going up and 
to the right) from its anchor point, and two other such rays (from other points) that limit it from 
the right and from the top. This immediately implies the existence of an exact algorithm for the 
optimization problem running in exponential time, based on brute force enumeration. The next 
lemma also shows that the greedy algorithm and brute force enumeration of permutations yields 
yet another exact algorithm for the optimization problem. 



Lemma 10 For every finite point set S C [0, 1]^ and every Pareto optimal anchored packing 
R = {r{s) : s £ S}, there is a permutation vr for which a greedy algorithm (with some tie breaking) 
computes R. 
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Proof. Let S C [0, 1]^ be a finite set, which may not contain the origin. Since every r{s) is Pareto 
optimal, it is a greedy choice assuming that s is the last point in the order vr. Suppose that s G S 
is the last point in a permutation vr. If the smaller problem S \ {s} with R \ {r{s)} is not Pareto 
optimal, then there is a point s' €z S for which we could choose a larger rectangle anchored at s', 
which intersects r(s) only. In this case, either s dominates s' or the ray shot from s' vertically 
up (resp., horizontally right) hits the lower (resp., left) side of rectangle r{s). This motivates the 
definition of a binary relation over S, which is an extension of the dominance order. Let s' ~< s 
if either s dominates s' or an axis-aligned ray shot from s' hits the lower or left side of rectangle 
r{s). It is not difficult to see that this is a partial order over S. If s G S is a minimal element in 
the poset {S, -<), then the rectangles in R\ {r{s)} are still Pareto optimal for the anchors S \ {s}. 
Now let TT be the reverse order of any linear extension of this partial order. □ 

We have shown (our main result) that for any set of n points in the unit square U = [0, 1]^, 
one can find a set of disjoint empty rectangles anchored at the given points and covering more 
than 9% of U. The same conclusion holds for points in any axis-aligned rectangle V instead of U, 
since it is straightforward to use an affine transformation to map the input into the unit square. 
Concerning the bound obtained, a sizable gap to the conjectured 50% remains. While certainly 
small adjustments in our proof can lead to improvements in the bound, obtaining substantial 
improvements probably requires new ideas. 
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